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Radar  Detection  of  Turbulence  in  Thunderstorms 


1.  INTRODUCTION 


The  detection  of  turbulence  in  thunderstorms  by  remote  methods  has  been  a 
subject  of  investigation  for  many  years.  Early  studies,  relying  on  observations 
from  penetrating  aircraft,  attempted  to  correlate  turbulence  encounters  with 
radar  storm  features.  The  general  conclusion  was  that  the  frequency  of  occur¬ 
rence  of  moderate  to  severe  turbulence  increased  markedly  with  increasing 
storm  radar  echo  intensity  (dBZ),  with  40  dBZ  a  rough  threshold  value  for 
development  of  such  turbulence  somewhere  in  the  storm.  Further  observation 
indicated  that  severe  turbulence  was  widely  distributed  in  storms  whenever  the 


intensity  indicated  damaging  hail  (about  50  dBZ),  with  turbulence  encounters 
outside  high  intensity  cores  being  as  frequent  as  those  within.  *  This  was  particu¬ 


larly  evident  within  thunderstorm  lines  where  severe  turbulence  was  frequently 
encountered  in  low  echo  intensity  regions  between  individual  storms.  The  con¬ 
clusion  was  that  storm  echo  intensity  was  an  indicator  of  the  presence  of  turbu¬ 
lence,  but  it  did  not  reliably  measure  its  strength  nor  locate  its  position. 

In  light  of  these  results  and  with  continued  development  in  Doppler  radar 
techniques,  emphasis  was  placed  on  finding  kinematic  signatures  of  hazardous 


(Received  for  publication  30  March  1981) 
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zones  within  thunderstorms.  Armstrong  and  Donaldson  developed  the  Plan 

Shear  Indicator  (PSI),  which  was  able  to  detect  regions  of  strong  shear  of  the 

3 

radial  wind.  Lee  and  Kraus,  by  comparing  PSI  and  penetrating  aircraft  ob¬ 
servations,  found  good  correlation  between  regions  of  strong  shear  of  the  radial 

4 

wind  and  severe  turbulence  events.  Radar  data  of  Frisch  and  Strauch  suggested 
that  high  turbulence  intensity  may  not  necessarily  be  correlated  with  zones  of  high 
shear  of  the  radial  wind,  but  when  present,  may  be  found  in  regions  exhibiting 
large  Doppler  spectrum  variance.  Later  Lee''  also  found  good  correlation  be¬ 
tween  regions  of  severe  turbulence  and  large  spectrum  variance.  Thus,  present 
attempts  are  being  directed  at  using  Doppler  spectrum  variance  as  a  turbulence 
indicator. 

The  turbulent  contribution  to  Doppler  spectrum  variance  is  a  function  of  the 
turbulence  intensity,  range  of  the  observation  volume  from  the  radar,  and  precip¬ 
itation  environment.  Turbulence  intensity  may  be  related  to  the  rate  per  unit  mass 
at  which  turbulent  kinetic  energy  is  transferred  from  larger  to  smaller  scales 
(eddy  dissipation  rate),  and  maximum  eddy  size  (turbulence  outer  scale).  The 
manner  in  which  all  these  effects  combine  is  not  well  known,  thus  the  interpreta¬ 
tion  of  this  indicator  is  not  yet  completely  understood.  The  aim  of  this  research 
is  to  clarify  the  relationship  between  Doppler  spectrum  variance  and  these  param¬ 
eters,  and  to  develop  a  more  reliable  radar  turbulence  detection  technique. 

This  report  focuses  on  results  of  theoretical  investigations  and  data  analyses 
performed  during  the  period  May,  1979  to  the  end  of  December  1980.  Primary 
topics  discussed  are:  1)  the  response  of  a  distribution  of  precipitation  particles 
to  various  ranges  of  scales  of  turbulent  motion;  2)  the  relationship  between  eddy 
dissipation  rate  and  Doppler  spectrum  variance  as  a  function  of  the  turbulence  scale 
regime  and  precipitation  environment;  3)  the  distribution  of  local  turbulence 
regimes  within  thunderstorms  as  delineated  by  aircraft  data;  and  4)  the  correla¬ 
tion  of  selected  aircraft  and  radar  data.  Finally,  since  one  of  the  ultimate  goals 
of  this  work  is  to  devise  a  methodology  for  turbulence  detection  by  incoherent 
radar,  a  possible  method  by  which  non-Doppler  radar  may  estimate  wind  shear 

2.  Armstrong,  G. ,  and  Donaldson  Jr. ,  R.  (1969)  Plan  shear  indicator  for  real¬ 

time  Doppler  radar  identification  of  hazardous  storm  winds,  J.  Appl. 

Meteor.  8:370-383. 

3.  Lee,  J.  T. ,  and  Kraus,  M.  (1975)  Plan  shear  indicator  and  aircraft  measure¬ 

ments  of  thunderstorm  turbulence:  experimental  results,  Preprints,  10th 
Radar  Meteorology  Conference,  Amer.  Meteor.  Soc.,  Boston,  pp.  337-340. 

4.  Frisch,  A.S.,  and  Strauch,  R.G.  (1970)  Doppler  radar  measurements  of 

eddy  dissipation  rates  in  a  northeastern  Colorado  convective  storm, 

J.  Appl.  Meteor.  15:1012-1017. 

5.  Lee,  J.T.  (1977)  Application  of 'Doppler  Weather  Radar  to  Turbulence 

Measurements  which  Affect  Aircraft.  Final  Report  FAA-RD-77-145  to 
Systems  Research  and  Development  Service,  FAA,  Washington,  D.  C. 
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2.  I  III.OKK1 1<  AL  l>\  ESTIMATIONS 


Pile  response  of  precipitation  particles  to  turbulent  air  motions  is  of  funda¬ 
mental  importance  since  the  Doppler  radar  velocity  is  generally  equated  to  the 
radial  component  (the  component  along  the  radar  viewing  direction)  of  the  envi¬ 
ronmental  wind.  However,  the  Doppler  power  spectrum  is  formed  from 
measurements  of  the  radial  velocities  of  the  precipitation  particles  within  the 
radar  pulse  volume,  with  the  spectral  power  at  each  spectral  velocity  value 
proportional  to  the  sum  of  the  sixth  powers  of  the  diameters  of  the  precipitation 
tracers  moving  with  that  radial  velocity  value.  It  is  suspected  that  a  sizable 
portion  of  Doppler  spectrum  power  may  be  supplied  by  tracers  which,  because 
of  their  large  inertia  (for  example,  large  raindrops),  do  not  follow  turbulent  air 
fluctuations  faithfully.  If  this  is  so,  then  an  incorrect  estimate  of  the  turbulent 
air  motion  (for  example,  Doppler  spectrum  mean  and  variance)  is  obtained. 

Thus,  one  is  interested  in  determining  the  response  of  a  distribution  of  particles 
to  various  scales  of  motion.  This  may  be  obtained  by  appropriately  summing  the 
individual  responses  of  the  particles  comprising  the  distribution. 

Consider  a  single  particle.  Its  equation  of  motion  may  be  written  as 


37tDp(VD  -  Va> 


CDNRe 

24 


(1) 


where  V^,  p^,  and  D  are  the  mass,  velocity,  density,  and  diameter, 

respectively,  of  the  particle.  Air  density  is  p  m^  is  the  mass  of  displaced  air, 
and  is  the  undisturbed  air  velocity  in  the  vicinity  of  the  particle.  Kinematic 
viscosity  is  p,  gravitational  acceleration  is  g,  and  Cp  and  (=2  Re)  are  the 
drag  coefficient  and  alternate  Reynolds  number  of  the  particle.  The  coefficient 
CdNr^/24  is  a  catchall  term  and  is  actually  related  to  fluid  acceleration  and 
past  particle  motion.  ®  This  equation  simply  relates  the  net  force  acting  on  the 
particle  to  the  pressure,  buoyancy-corrected  gravity,  and  drag  terms. 

Equation  (1)  may  be  written  as 


6.  Pruppacher,  H.  R.  and  Klett,  J.D.  (1978)  Microphysics  of  Clouds  and 
Precipitation,  D.  Reidel  Publishing  Co. ,  Boston. 
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where 
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(3) 


Actual  solution  of  Eq.  (2)  requires  application  of  numerical  methods,  since 
Cj-j  is  dependent  on  (V^  -  V  ).  To  obtain  a  general  analytic  solution  for  the  par¬ 
ticle  response  it  is  necessary  to  make  a  simplifying  assumption  concerning  the 
drag  force.  For  this  purpose  it  is  assumed  that,  on  the  average,  the  particle  is 
settling  at  its  terminal  fallspced  (V7^,).  The  drag  force  in  Eq.  (2)  may  then  be 
replaced  bv  its  average  value. 

7 

Following  Stackpole,  the  following  relation  is  then  applicable: 


Neglecting  particle  growth  and  the  variation  of  air  density  with  height,  r  is  a  con¬ 
stant  for  each  particle.  Note  that  this  assumption  means  that  the  drag  coefficient 
is  replaced  by  the  value  attained  when  the  particle  is  settling  at  its  normal  termi¬ 
nal  fallspeed  in  still  air,  and  thus  it  is  a  slightly  averaged  particle  response  that 
may  be  estimated. 

In  applying  Eq.  (2)  to  a  turbulent  field,  one  generally  considers  turbulence  as 
composed  of  a  random  collection  of  eddies  having  a  range  of  scales  (for  example, 
wavelengths  X).  Mathematically,  one  may  view  this  collection  as  a  series  of 
harmonic  functions  of  varying  radian  frequency  (u)  or  wavenumber  (|k]  -  2tt/X), 
having  random  amplitude  and  combining  to  produce  the  observed  turbulent  struc¬ 
ture.  The  response  of  a  particle  to  a  turbulent  gust  may  therefore  be  considered 
the  result  of  a  superposition  of  the  particle's  responses  to  this  collection  of  har¬ 
monic  functions,  or  eddies. 

The  form  of  the  par  ticle  response  to  a  single  such  function  may  be  seen  by 
considering  one  dimension  (horizontal)  and  introducing  an  air  velocity 


7.  Stackpole,  J.  D.  (lPfil)  The  effectiveness  of  raindrops  as  turbulence  sensors, 

Proceedings  Pth  Weather  Radar  Conference.  Amer.  Meteor.  Soc. ,  Boston,  ' 

pp.  212-217. 
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inti.  Eq.  (2).  The  particle  velocity  is  found  to  In 


VD  Vo 


f  ■>  ,  ■>  i  v  i  A  i(u-T+ATANuT(, >  -1» 

Ul  +  w“T-p“r  ~/(l  +  u.“ T ”)  -  J e 


A  e 


■T  it 


CO 


where 
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Note  that  the  solution  for  the  vertical  dimension  has  an  additional  vertical  accel¬ 
eration  term  g.  The  first  term  shows  that  the  particle  velocity  is  reduced  in 
amplitude  and  phase  shifted  relative  to  the  air  velocity.  The  second  term  is  a 
transient  term  and  describes  the  time  required  to  achieve  this  reduced  amplitude, 
phase  lagged,  state  of  motion.  The  phase  shift  alone  does  not  represent  any 
ener  gy  loss  for  the  particles.  However,  if  the  particles  are  entering  a  region  of 
greater  turbulent  air  energy  (for  example,  a  r  egion  dominated  by  larger  eddies), 
then  the  transient  period  represents  a  time  when  the  particles  have  less  energy 
than  they  will  exhibit  once  they  have  achieved  their  new  state  of  oscillatory  motion. 
On  the  other  hand,  if  the  particles  are  entering  a  region  of  decreased  turbulent 
air  energy  dominated  by  smaller  scale  eddies,  then  the  particles  have  excess 
energy  during  the  transient  period.  Thus,  on  the  average,  the  transient  effects 
may  not  represent  any  net  turbulence  energy  loss  when  compared  to  the  hypothetical 
case  where  the  particles  have  no  transient  period  at  all  (where  they  attain  their 
new  state  of  motion  instantaneously). 

Furthermore,  assuming  the  particles  to  be  raindrops  having  a  range  of 
terminal  fallspeeds  of  1  m/sec  to  9  m/sec,  the  transient  term  in  Eq.  ((5)  indicates 
that  it  would  take  0.  25  sec  to  2  sec,  or  0.  25  m  to  18  m  distance  fallen,  to  pass 
through  the  transient  period.  These  lengths  represent  only  the  very  lowest  energy 
portion  of  the  turbulence  regimes  of  interest  here  (outer  scale  lengths  of  ^500  m). 
Thus,  the  net  transient  energy  loss  should  be  very  small  indeed.  Therefore, 
since  the  mathematical  techniques  involved  here  deal  only  with  the  average  tur¬ 
bulent  motion  of  a  distribution  of  particles  in  relation  to  the  average  turbulent 
air  energy,  only  the  amplitude  response  term  in  Eq.  (6)  will  be  considered. 

It  is  necessary  to  map  this  response  function  from  dependence  upon  frequency 
to  turbulence  scale,  since  frequency  is  dependent  upon  the  relative  speed  of  the 
atmospheric  structure  and  the  sensing  instrument  (here  the  particles),  and 
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therefore  is  an  artificial  turbulence  parameter.  In  a  manner  similar  to  Taylor's 
hypothesis,  one  may  form  a  relationship  between  turbulence  wavenumber,  r  adian 
frequency,  and  particle  terminal  fallspeed  as 

u  :  VT  •  K  =  VtKz  ,  (7) 

resulting  in  the  response  function  taking  the  form 

VD/Va  1  l/d  +  K-V^./g2)1'2  ,  (B) 

where  o  ~  10  2  and  uitp  <  1  have  been  assumed.  This  last  condition  is  met  for 

o  o 

all  rain  and  snow  environments  except  for  those  combining  large  particle  size 
and  particle  observed  scales  less  than  about  0.0a  m.  The  energy  content  of  the 
turbulent  field  at  these  scales  is  very  small  in  relation  to  the  larger  scales  of 
interest  here  and  need  not  be  considered  further, 

7 

This  relationship  was  first  obtained  bv  Stackpole  and  states  that  as  b  • 
particle  size  increases  (V^.  increases)  the  response  to  turbu’>  •->  motions  having 
any  given  apparent  scale  A  (=  Kz/2tr)  decreases.  Conversely,  a  particle  of  given 
size  (V^,  constant)  exhibits  decreasing  response  to  decreasing  scale.  The  term 
"apparent  scale"  has  been  used  to  signify  that  the  particles  are  responding  to  the 
field  which  is  the  projection  of  the  three  dimensional  turbulent  field  upon  the 
z-axis  and  that  there  may  not  be  a  one  to  one  correspondence  of  features  in  the 
three  dimensional  field  to  those  observed  in  the  one  dimensional  field. 

To  assess  the  effect  imperfect  particle  response  has  upon  the  measurable 
radar  parameters,  it  is  necessary  to  modify  the  basic  pertinent  radar  relations. 
In  all  that  follows,  the  wind  environment  is  considered  fully  turbulent,  isotropic, 
and  homogeneous.  It  is  further  assumed  that  the  effects  of  Doppler  spectrum 
broadening  factors  are  additive  and  thus  each  contribution  may  be  viewed 
separately.  Here  only  the  purely  turbulent  contribution  will  be  investigated. 

In  general  form  the  Doppler  spec  trum  variance  may  be  written  as 


2  —  o 

VAR  -  V  -  V“ 
r>  o 


where  V  is  particle  radial  velocity  and  the  overbar  indicates  averaging  >.vi  i  ');<■ 

p 

Doppler  power  spectrum.  This  relation  states  that  the  spectrum  variance  equai. 
the  averaged  square  of  particle  radial  vcloi  ity  less  the  squar-  if  -he  averaged 
particle  radial  More  spe-iftcelb-,  'res  equaMon  itf»v  be  written  as 
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JJ V». 

1  *  *  p 


(R)q(Vp,  R)l(r,  K)dV  dK 


VA  H 


cjf  q(Vp,  iT)I(r%  t?)dVpdR 


C j//  Vp(R)q(Vp,  R)l{r,  R)dVpdR 
R _ 

cjf  q(Vp,  R)I(r*  R)dVpdR 
R 


(10) 


showing  more  clearly  the  backscattered  power  averaging  and  pulse  volume  aver¬ 
aging  over  available  particle  radial  velocities.  Here  is  the  radar  constant, 
and  q(Vp,  R)  is  the  reflectivity  density  of  the  particles  in  the  volume  element  dR 
(related  to  the  sixth  power  of  particle  diameter,  or  melted  diameter  for  snow) 
located  R  from  the  pulse  volume  center  located  at  r  (relative  to  the  radar),  and 

having  radial  velocities  ranging  from  V  to  V  +  dV  .  The  term 

P  P  P 

C',q(V  ,  R)I(r,  R)dRdV  is  the  power  backscattered  by  these  particles  to  the  radar, 
and  I(r,  fi)  is  the  two-way  beam  illumination  pattern  which  is  assumed  normalized. 
The  integral  over  R  is  the  integral  over  the  pulse  volume,  and  the  Vp  integral  is 
over  the  total  unambiguous  velocity  range  of  the  radar. 

Generally,  the  precipitation  backscattered  power  varies  with  position  within 
the  pulse  volume;  however,  to  incorporate  the  effect  of  imperfect  response  and 
make  the  above  relation  tractable,  a  simplifying  assumption  is  required.  It  is 
assumed  that  during  the  time  the  radar  samples  a  given  region  to  form  the  sam¬ 
ple  data  set  (used  to  form  the  Doppler  power  spectrum),  a  characteristic  dis¬ 
tribution 

N(D)  =  N  e"AD  (11) 

o 

of  particles  is  swept  past  each  point  within  the  pulse  volume.  Here  N(D)  dD  is 
the  number  of  particles  per  unit  volume  in  the  size  range  D  to  D  +  dD.  In  essence, 
this  states  that  the  reflectivity  is  constant  throughout  the  pulse  volume. 

The  air  velocity  at  any  given  point  within  the  pulse  volume  is  the  result  of  a 
superposition  of  various  sized  eddies  to  which  each  particle  in  the  distribution  has 
a  definite  response  determined  by  its  size.  Thus,  it  is  assumed  that  for  a  given 
eddy  distribution,  each  point  within  the  pulse  volume  has  a  well  defined 
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relationship  between  particle  Doppler  speed  and  size.  This  allows  the  dependence 
upon  particle  velocity  to  be  changed  to  dependence  upon  particle  diameter 

q(Vp,  R)dVp  -»  q(D)  dD  (12) 

and  allows  Eq.  (10)  to  be  rewritten  as 


VAR 


/  Vpd(^)I(,r-  ^)d[*  ' 
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2 

/  Vpd(K>I(»V  KMR 

R 


where 


(13) 


Jv2d(  R)q(D)dD 
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q(D)dD 


and 
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J  V^(R)q(D)dD 
j  q(D)dD 


(15) 


where  V^d?)  is  the  radial  velocity  of  a  given  drop  of  diameter  D  located  at  the 
position  R  in  the  pulse  volume. 

This  is  essentially  the  classical  definition  of  Doppler  spectrum  variance. 
However,  the  significant  difference  here  is  that  V  ^(R)  is  the  reflectivity-weighted 
imperfect  response  of  the  particles,  and  thus  incorporates  the  concept  that  during 
the  sampling  period,  at  a  given  point  in  space,  the  various  size  particles  compos¬ 
ing  the  particle  distribution  exhibit  different  velocities  from  each  other-  and  the 
single-valued  air  velocity  at  that  point.  With  an  ensemble  (denoted  by  (  ' )  of 
observations,  and  noting  that  ensemble  averaging  mathematically  is  a  summation 
and  -an  be  brought  within  the  integral,  Eq.  (13)  becomes 

<VAR)  =  (lV2pd))  -<(Vpdl2N  <V2dN  -VAR([Vpd!)  ,  <1B) 

I 

1 


where  now  ihe  bracket  indicates;  averaging  over  the  pulse  volume  (note  that 
<[ Vpd]  > .  the  average  Doppler  spectrum  mean  velocity  in  the  turbulent  field,  is 
zero).  This  states  that  the  average  Doppler  spectrum  variance  equals  the  aver¬ 
age  one-dimensional  reflectivity-weighted  particle  turbulent  motion  variance 
minus  the  average  variance  of  the  Doppler  mean  velocity.  This  is  a  modification 

g 

of  the  relation  of  Rogers  and  Tripp.  It  is  important  to  note  here  that  this  radar 
mean  particle  turbulent  variance  is  not  the  air  turbulent  variance,  but  the  reflec¬ 
tivity-weighted  response  of  the  particle  distribution  to  the  air  turbulent  motions. 
Similarly,  the  Doppler  mean  velocity  term  also  includes  the  imperfect  response 
function. 

To  incorporate  the  response  functions  it  is  useful  to  transform  Eq.  (lfi)  into 

Q 

wavevector  space.  Srivastava  and  Atlas"'  have  shown  that  the  average  variance 
of  the  Doppler  spectrum  (Eq.  (l(i|  ,  for  the  case  of  a  uniform  distribution  of 
perfect  tracers)  is  given  by 


(VAR)  J  <t,u(K)(l  -  ^(K))  dK 


where  Kis  the  ihree-dimensional  turbulent  wavevector,  <J^(K)  -  (2»)l,Fj(K)F*(K), 
where  Fj(K)  is  the  Fourier  transform  of  Hr,  R),  and  c>jj(K)is  the  turbulent  power 
density  function  of  the  radial  (1)  air  velocity  component.  Equation  (17),  however, 
does  not  include  the  particle  response  functions.  Referring  back  to  Kqs.  (13) 
through  (15),  the  new  modified  form  of  this  relation  becomes 


(VAR)  -  J  <Ju(K)(R,)(Kz)  -  RjdOOjfK))  dK 


CJ  /  "'A,U  D"  f  2  dD 

Z  •'  (1  +  K^(aDb)4/g2)  j  J 
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9.  Srivastava,  R.C.,  and  Atlas,  D.  (1974)  Effect  of  finite  radar  pulse  volume  on 
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"AD  I)'’  dD 


The  particle  fallspeed  has  been  replaced  by 
VT  -  a  Db  , 


(20) 


(21) 


where  D  is  particle  diameter  (or  melted  particle  diameter)  and  a,  b  are  constants. 
The  reflectivity  weighted  (D*’  weighting)  imperfect  response  function  is  now 
explicitly  shown.  The  terms  R.,(K  ),  R,(K  )  are  the  particle  turbulent  motion 
variance  and  turbulent  velocity  weighting  factors,  respectively. 

In  Eq.  (18),  ^(K)  essentially  describes  the  actual  one-dimensional  ’urbulent 

air  motion  intensity  as  a  function  of  turbulent  scale.  Multiplication  of  tins  by 
R^(K^)  results  in  the  description  of  maximum  radar  measurable  one-dimensional 
particle  motion  varianc  e  as  a  function  of  turbulent  scale.  Multiplying  o^fK)  by 
Rl(Kz)6j(K)  describes  how  much  of  this  radar  measurable  precipitation  motion  is 
mapped  into  fluctuation  of  the  mean  Doppler  veloc  ity.  In  both  cases,  the  action 
of  the  response  function  is  to  pass  along,  relatively  unscathed,  the  largo  scale 
turbulent  motions  to  the  precipitation,  while  severely  limiting  the  small  scale 
turbulent  motions.  The  average  Doppler  spectrum  variance  is  simply  the  differ¬ 
ence  between  the  maximum  one-dimensional  radar  measurable  precipitation  motion 
variance  and  the  average  variance  of  the  fluctuation  of  the  Doppler  mean  velocity 
about  the  ensemble  average  of  zero.  This  is  the  desired  result. 

There  are  a  few  points  to  note  tn  Eq.  (18).  First,  the  results  are  dependent 
only  on  the  slope  (A)  of  the  particle  distribution  function  and  not  on  the  actual 
particle  .  oncentration  number  X  .  Secondly,  if  the  pulse  volume  is  decreased 
in  size  to  a  point  (in  which  case  the  beam  filter  function  djr,  R]  behaves  as  a 
delta  function,  for  which  bj[K|  becomes  unity  over  all  ft,  ),  there  exists  a  min¬ 
imum  variance  which  must  always  be  observed.  Its  value  is  given  by 


<VAlDmin  -  /  <*>u(K)(R2(Kz)  -  R^K.,))  dK  .  (21) 

K 

Numerical  calculations  show  that  this  term  is  small  and  reflects  the  fact  that 
the  two  weighting  functions  are  nearly  identical  in  net  effect.  The  existence  of 
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this  minimum  variance  is  expected,  since  the  point  in  space  (nost  representing  the 
pulse  volume)  is  observed  for  a  finite  am<  unt  of  time,  during  which  a  distribution 
of  particles  having  differing  radial  velocities  (because  of  varying  response  charac¬ 
teristics)  sweeps  bv  and  contributes  to  the  Doppler  power  spectrum. 

The  variation  of  these  response  functions  with  wavenumber  component  K  is 
shown  in  Figures  la  and  lb  for  var  ious  particle  distr  ibution  slopes.  The  particles 
are  assumed  to  be  raindrops  with  the  constants  in  Kq.  (21)  set  at 

a  -  1(190  cm  *'  ^  and  b  ---  0.  (1 

These  figures  show  that  the  reflectivity-weighted  precipitation  response  is  good 
for  turbulent  motions  having  an  apparent  scale  (measured  along  the  '/-direction) 
larger  than  about  250  m,  with  virtually  no  response  to  scales  of  motion  smaller 
than  about  10  m.  As  expected,  better  response  is  obtained  in  light  rain 
(A  -  ‘>0  cm  1 )  than  in  heavy  rain  (A  =  20  cm  *)  environments.  It  should  also  be 
remembered  that  the  response,  although  excellent,  is  not  perfect  for  scales  greater 
than  1  km. 

To  foster  a  better  appreciation  for  the  effects  of  these  i  esponso  fun<  tiotis,  it 
is  useful  to  compare  the  maximum  three-dimensional  precipitation  motion  variance 
measurable  by  radar  ((VAR^  )  and  that  actually  contained  by  the  particles 
((VAR''rn).  The  radar  measurable  quantity  demonstrates  R.^iK^)  and  could  be  the 
sum  of  the  mean  Doppler  spectrum  variance  and  average  variance  of  the  mean 
Doppler  velocity  from  three  radars  observing  the  same  location  from  orthogonal 
directions.  The  mechanical  precipitation  variance  demonst rates  Kq.  (H).  The 
maximum  three-dimensional  mechanical  turbulent  motion  variance  is  given  by 

<VAR)m  =  <u2  +  v2  +  w2)m  =  J  E(K)W(K)  dK  ,  (23) 

where  the  air  energy  density  E(K)  is  assumed  given  by  the  Kolmogorov  relation 

E(K)  =  C  e2/,3K'5/3  (24) 

in  which  C  is  a  known  constant,  e  is  the  eddy  dissipation  rate,  and  K  =  |k|.  The 
particle  distribution  response  is 


W(K)  =  2 J  [[gN(D)atan(KV2/g)]/KV2  j  dD^  f  N(D)  dD 


(25) 
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this  minimum  variance  is  expected,  since  the  point  in  space  (nov\  representing  the 
pulse  volume)  is  observed  for  a  finite  amount  of  time,  during  which  a  distribution 
of  particles  having  differing  radial  velocities  (because  of  varying  response  charac¬ 
teristics)  sweeps  by  and  contributes  to  the  Doppler  power  spectrum. 

The  variation  of  these  response  functions  with  wavenumber  component  is 
shown  in  Figures  la  and  lb  for  various  particle  distribution  slopes.  The  partirles 
are  assumed  to  be  ra indrops  w  ith  the  constants  in  Kq.  (21)  set  at 

a  -  1  f> 90  cm^‘4  and  b  =  0.  6 

These  figures  show  that  the  reflectivity-weighted  precipitation  response  is  good 
for  turbulent  motions  having  an  apparent  scale  X  ^  (measured  along  the  /.-direction) 
larger  than  about  250  m,  with  virtually  no  response  to  scales  of  motion  smaller 
than  about  10  m.  .As  expected,  better  response  is  obtained  in  light  rain 
(A  -  ''0  cm  *)  than  in  heavy  rain  (A  “  20  cm  *)  environments.  It  should  also  be 
remembered  that  the  response,  although  excellent,  is  not  perfect  for  scales  greater 
than  1  km. 

To  foster  a  better  appreciation  for  the  ••ffects  of  these  response  functions,  it 
is  useful  to  compare  the  maximum  three-dimensional  precipitation  motion  variance 
measurable  by  radar  ((VAFF  )  and  that  actually  contained  bv  the  particles 
((VARSm).  The  radar  measurable  quantity  demonstrates  )  and  could  be  the 

sum  of  the  mean  Doppler  spectrum  variance  and  average  variance  of  the  mean 
Doppler  velocity  from  three  radars  observing  the  same  location  from  orthogonal 
directions.  The  mechanical  precipitation  variance  demonstrates  Eq.  (8).  The 
maximum  three-dimensional  mechanical  turbulent  motion  variance  is  given  by 

<vAR)m  3  <ti2  +  v2  +  w2>m  =  J  E(K)W(K)  dK  ,  (23) 

where  the  air  energy  density  E(K)  is  assumed  given  by  the  Kolmogorov  relation 

E(K)  =  C  e2/,3K‘5/3  (24) 

in  which  C  is  a  known  constant,  e  is  the  eddy  dissipation  rate,  and  K  =  |k|.  The 
particle  distribution  response  is 


W(K)  =  2j  {  [gN(D)atan(KV^,/ g)j  /KV^,}  dD  j  J  N(D)  dD 


(25) 
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where  the  particle  number  density  and  fallspeed  are  given  by  Eqs.  (11)  and  (21), 
respectively.  Figure  2  shows  this  variance  (VAR)  ,  normalized  to  the  ait- 
motion  variance  (VAR)^  where 


(VAR)a  =  2  J  E(K)  dK 


(26) 


0  4  - 1 

The  constant  parameters  in  Eqs.  (23)  through  (26)  are  a  =  1690  cm  '  sec 

b  =  0,  6  D  0.  005  cm,  D  =  0.  5  cm,  and  K  =  2ircm  *  (corresponding 

nun  ’  max  max  b 

to  a  minimum  turbulence  scale  of  1  cm).  Parameters  K  (~  2tr/A^  )  correspond¬ 
ing  to  a  turbulence  outer  scale  A  and  the  particle  distribution  slope  A  are  the  inde¬ 
pendent  variables. 

This  figure  shows  that  the  mean  turbulent  variance  of  the  particles  is  typically 
95  percent  (or  greater)  of  the  air  value,  for  all  precipitation  environments  (for 
example,  heavy  rain  A  =  10  cm  ^  to  light  rain  A  =  60  cm  *)  and  outer  scales  (A^) 
greater  than  500  m. 

This  figure  demonstrates  that  eve:  .o  .gh  individual  particles  (for  example, 
large  drops)  may  net  respond  well  in  smal.  scales  of  motion,  the  actual  mean 
turbulent  velocity  variance  of  the  particle  distr  ibution  is  generally  a  good  measure 
of  the  atmospheric  value.  Unfortunatr  tv,  the  radar  does  not  weight  the  particles 
identically,  as  this  calculation  assumes,  but  rather  weights  the  particles  according 
to  the  sixth  power  of  drop  diameter. 

The  maximum  three-dimensional  radar  measurable  turbulent  variance  is 
given  by 

(VAR)  =(u2  +  v2  +  w2)  =  f  E(K)W,(K>  dK  ,  (27) 

x  rm  rm  J  1 

where 


W  j(K)  =  2  J,{[gN(D)D6atan(KV2/g)]  /KV^,}dD  j  J  N(D)Df 


dD 


(28) 


This  quantity,  normalized  to  the  actual  turbulent  air  variance,  is  shown  in 
Figure  3.  It  is  observed  that  if  the  outer  scale  (X  )  extends  to  a  kilometer  or  more, 
then  the  radar  can  measure  90  to  97  percent  of  the  air  motion  variance  in  all  pre¬ 
cipitation  environments.  For  shorter  outer  scales,  however,  particularly  in  heavy 
rain  situations,  the  radar  estimate  may  be  in  significant  error.  As  an  example, 
for  A  =  10  cm  1  and  an  outer  scale  of  250  m,  the  total  radar  measurable  variance 
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is  only  about  72  percent  of  the  air  value.  The  striking  difference  between  Figures  2 
and  3  results  from  the  strong  power  contribution  to  the  Doppler  spectrum  from  large 
particles  which  exhibit  poor  response  character ist ics,  and  demonstrates  that  the 
usual  assumption  of  precipitation  being  perfect  tracers  of  the  air  motion  can  lead 
to  radar  turbulence  energy  estimates  having  significant  error.  The  results  demon¬ 
strate  that  imperfect  particle  response  and  the  finite  outer  scale  of  turbulence  must 
be  taken  into  account  if  reasonable  estimates  of  turbulence  parameters  are  to  be 
obtained  by  radar  methods. 

It  was  stated  that  the  available  radar  measurable  precipitation  motion  variance 
is  partitioned  between  the  average  Doppler  spectrum  variance  and  the  average  var¬ 
iance  of  the  fluctuating  Doppler  velocity.  This  effect  may  best  be  illustrated  by 
graphical  presentation  of  the  two  terms  forming  the  right  side  of  I.'q.  (18).  For 
mathematical  ease,  the  viewing  direction  is  chosen  as  vertical.  The  radial  direc¬ 
tion  (1)  now  becomes  the  z  direction.  The  resulting  average  Doppler  spectrum  var¬ 
iance  relation  (Eq.  [18J )  may  be  written  as 


<VAR>  J  <*zz<Kz>R2(Kz> 


*zzf(Kx)Hl<Kz»  dKz 


(29) 


where 


*zz<Kz>  II  <*,zz^d*Scd^y  ^0) 

♦  zzW  =  ff*  zz(l?)<¥^)dKxdKy  (31) 


and 


*zz<K) 


(32) 


Here  <J>  z(Kz)  's  the  true  turbulent  air  one -dimensional  longitudinal  power  density 
spectrum  function  along  the  z-direction.  Similarly,  ^zzf(Kz)  represents  the  spec¬ 
trum  that  would  result  from  power  spectrum  analysis  of  Doppler  radial  velocities 
obtained  from  a  set  of  successive  range  bins  (close  pnough  so  that  the  pulse  vol¬ 
ume  dimensions  may  be  considered  constant)  along  the  radial  at  some  range  r 
from  the  radar  in  an  environment  of  a  uniform  distribution  of  perfect  tracers. 
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Tile  ■  'lit it  I'uta'c  t!  .heled  pt)  is  the  true  turliuletlt  all  -.pi  a’jun,  o  (K.  )  [Fq.  ('iO)l 
.  id  ihcoi -o- irally  i  epiv:  ••■its  the  sport  run:  'hat  would  !■••  .blam.-d  from  power  spo.  • 
t  rum  analysis  of  the  Doppli  r  velocities  from  sifees  ••:%*»•  p.  ar.t  (pt)  sized  pulse  vol¬ 
umes  along  the  z-diree’itiR  in  an  envituniien'  of  pe.fee-  traiers  (an).  The  low.-t 
i  n:  vis  (labeled  a,  20,  40,  00,  120)  are  tin-  i  orrespondmg  radar  pulse  volume  fil  — 
•ered  turbulent  air  (perfect  tracer)  spectra  <? V.,P(K  )  1 1  Iq .  (41)]  for  tlie  five  ranges 
mentioned,  The  curves  labeled  pt  A  10  iitr*  or  pt  A  •  > 0  cm  ^  represen'  the 
modifications  to  the  true  turbulent  air  spectra  for  rain  environments  bavins 
A  10  cm  ^  and  A  t >0  cm  *  and  represent  the  o  (K  MI./K  )  term  in  i.q.  (2!'). 

The  curve  labeled  20  km  A  10  cm  is  'he  turbulence  putter  spectrum  that  ..uld 
result  when  the  pulse  volume  dimensions  are  'hose  of  a  real  pulse  volume  located 
20  km  from  the  radar,  and  is  representative  of  the  true  radar  measurement.  Die 
imperfect  tracer  curves  for  ranges  greater  than  20  km  are  nearly  identical  to  the 
respective  perfect  tracer  (air)  curves  and  are  riot  drawn.  I  his  last  effect  simply 
reflects  the  better  response  of  the  precipitation  particles  to  the  dominant  contrib¬ 
uting  scales,  which  arc  increasing  with  increasing  pulse  volume  size.  These 
curves  show  the  reduction  of  turbulence  energy  at  small  apparent  spectral  scales 
(large  Kj)  due  to  imperfect  precipitation  response.  This  indicates  that  as  research¬ 
ers  reduce  the  observation  range  to  reduce  the  effects  of  pulse  volume  filtering, 
they  must  pay  greater  attention  to  the  effect  imperfect  precipitation  response  w  ill 
have  on  their  derived  turbulence  power  spectra. 

The  curves  also  show  that  the  "knee"  of  the  curves  (Figure  4b)  are  a  good 
indication  of  the  turbulence  outer  scale  length  for  this  isotropic  field  |Fq.  (32)1.  We 
will  next  return  to  the  task  of  extracting  useful  turbulence  information  from  the 
Doppler  spectrum  variance. 

The  Doppler  spectrum  variance  [Eq,  (29)1  is  the  areal  difference  between  the 
appropriate  point  and  range  spectral  curves  (Figure  4b  but  plotted  on  linear  scales). 
For  example,  in  a  perfect  tracer  (air)  environment  having  a  turbulence  outer  scale 
A  =1  km,  the  Doppler  spectrum  variance  for  the  case  where  the  radar  is  observ¬ 
ing  a  region  20  km  away  is  the  areal  difference  between  the  air  point  curve  and  air 
20  km  curve  in  Figure  4b.  Similarly,  in  a  rain  environment  having  A  r  10  cm  , 
the  Doppler  spectrum  variance  is  the  areal  difference  between  the  pt  A  -  10  cm  1 
curve  and  the  20  km  A  "  10  cm  *  curve.  It  should  be  noted  that  for  ranges  greater 
than  about  20  km,  the  pulse  volume  filtered  spectra  for  air  (#zzf)  may  be  substituted 
for  the  precipitation  environment  curve  with  negligible  error. 

To  determine  whether  the  interpretations  from  the  vertical  observations  can  be 
extended  to  the  horizontal  direction,  calculations  were  performed  to  investigate  the 
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partitioning  of  the  total  radar  measurable  motion  variance  ((u^s  +  (v^)  +  {w^') 
among  the  three  orthogonal  components.  The  analysis  shows  that  in  light  rain 
situations,  this  quantity  was  partitioned  equally  between  the  x  v  z  components. 

In  heavy  rain  environments,  however,  the  horizontal  variance  components  were 
noticeably  smaller  than  the  vertical  component  for  turbulence  outer'  scale  lengths 

2  9 

less  than  about  0.  2a  km  (for  example,  fu  )  ,( w“)  -  0.  9,r),  0.  91,  0.  88  for-  A  =  0.  5, 

'  ■  '  o 

0.2a,  0.  197  km  respectively).  This  is  a  real  effect  and  results  from  the  horizontal 
and  vertical  power  density  spectra  having  a  different  K  dependence,  while  the 
response  function  of  the  particle  disti  ibution  is  dependent  only  upon  K,  .  A  more 
exact  form  of  the  response  function  allowing  for  some  slight  dependence  upon  K„ 

K  would  presumably  retain  greater  equality  among  (he  three  eomponent  motion 
variances  to  shorter  outer  scale  values. 

To  determine  the  relevance  of  these  numbers,  consider  the  longitudinal 
power  density  spectra  that  would  he  obtained  by  looking  along  the  x  ur.d  / 
directions  in  a  turbulent  environment  with  A  i  km.  for  the  ease  of  pc-rfer* 
t  raeers  (air),  the  spectra  o  (K  )  and  o  (K  )  are  identical  to  the  air  point  ■  urve  in 
Figure  4b.  No'e  that  the  area  below  each  curve  would  tie  (u“  =  (w  .  Now, 

introduction  of  precipitation  will  produce  the  modified  point  curves  )H.,(Kz) 

and  ^  (K  )Rg(K  ).  If  A  10  cm  \  the  vertical  direction  curve  is  identical  to  the 

pt  A  -  10  cm"'  curve  in  Figure  4b.  However,  because  of  the  different  K,  depend - 

- 1  z 

ence  of  the  horizontal  turbulence'  spectrum,  the  pt  A  10  cm  curve  f. the 
horizontal  direction  will  be  slightly  different.  Similarly,  for  finite  observation 
ranges,  the  u  and  w  (Figure  4b)  turves  will  continue  to  exhibit  slight  differences. 
However,  the  near  equality  ((u  x  ’  (w  ))  for  cases  where  A^  >  0.2"  km  suggests 
that  the  differences  in  the  u  and  w  precipitation  environment  longitudinal 
power  spectral  curves  will  be  at  very  low  apparent  scale  (large  K^,  Kz)  values. 

For  finite  ranges,  pulse  volume  filtering  dominates  over  precipitation  response 
effects,  thus  the  difference  in  Doppler  spectrum  variances  results  primarily  from 
the  difference  in  the  modification  of  the  point  curves  which  is  known  to  be  small. 
These  results  indicate  that  solution  of  Eq.  (25)  yields  the  Doppler  spectrum  var¬ 
iance  that  would  be  measured  along  any  viewing  direction  in  any  precipitation  envi¬ 
ronment  exhibiting  turbulence  outer  scale  values  of  0.  25  km  or  greater. 

Numerical  evaluation  of  Eq.  (29)  was  performed  to  estimate  the  Doppler  spec¬ 
trum  variance  that  would  be  obtained  in  a  variety  of  precipitation  and  turbulent 
environments.  Figures  5a  and  5b  show  the  variance  (normalized  to  C  versus 

the  range,  and  the  ratio  of  pulse  volume  width  to  depth.  The  solid,  long  dash,  and 
short  dash  curves  represent  the  air  (perfect  tracer)  and  imperfect  tracer  (A  =  60, 

10  cm  *)  environments  respectively.  The  curve  marked  \q  =  oo  is  essentially 
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VAR 


Figure  5a.  Doppler  Spectrum  Variance  (VAR)  Normalized 
to  C  e2/3  Versus  Radar  Range,  Turbulence  Outer  Scale 
Lengths  are  =  1,  4,  oo  km.  Full  half-power  beamwidth 
is  1”  and  pulse  volume  depth  is  200  m.  Ratio  of  pulse  vol¬ 
ume  width  to  depth  ratio  at  top 
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Figure  5b.  Doppler  Spectrum  Variance  (VAR>  Normalized 
to  C  e  ■  J  Versus  Radar  Range,  Turbulence  Outer  Scale 
Lengths  are  A.Q  =  0.  5,  2,  8  km.  Full  half-power  beamwidth 
is  1°  and  pulse  volume  depth  is  200  m.  Ratio  of  pulse  vol¬ 
ume  width  to  depth  ratio  at  top 


Figure  4a  and  is  well  represented  by  a  corrected  form  of  the  analytical  approxima- 
tt.m  of  Frisch  and  C  lifford.  ^  It  is  noted  that  the  spectrum  variance  is  essentially 
independent  of  outer  scale  length  for  ranges  less  than  about  20  km.  This  simply  re¬ 
flects  that  at  short  ranges  the  various  turbulent  fields  are  nearly  indistinguishable 
to  the  relatively  small  (compared  to  outer  scale  size)  pulse  volume,  and  the  variance 
contribution  can  he  reasonably  estimated  from  the  corrected  form  of  the  approxima¬ 
tion  of  Frisch  and  C  lifford.  Significant  deviation  from  this  relation  occurs  when 
the  maximum  pulse  volume  dimension  approaches  the  outer  scale  value.  When  the 
maximum  pulse  volume  dimension  becomes  larger  than  the-  turbulence  outer  scale, 
nearly  all  turbulent  motion  is  mapped  into  Doppler  spectrum  variance.  As  radar 
range  increases  beyond  this  point,  the  Doppler  variance  becomes  constant,  inde¬ 
pendent  of  range. 

The  variation  due  to  precipitation  environment  is  significantly  smaller,  with 
increasing  dependence  on  rain  environment  occurring  with  decreasing  ou'er  scale 
size  and  decreasing  range.  These  features  can  more  easily  be  seen  bv  inverting 
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10.  Frisch,  A.S.  ,  and  C  lifford,  S.  F.  (1074)  A  study  of  convection  capped  by  a 
stable  layer  using  Doppler  radar  and  ncoustii  echo  sounders,  J.  Atmos. 
Sri.  3  1:  1 022-1  128. 
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Figure  6a.  Eddy  Dissipation  Rate  Normalized  to 
((VAR) /C)3/2  Versus  Radar  Range,  Turbulence 
Outer  Scale  Length  are  X.p  =  1,  oo  km.  Radar 
full  half-power  beamwidtn  is  1°  and  pulse  volume 
depth  is  200m.  Ratio  of  pulse  volume  to  depth  at 
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Figure  fib.  Eddy  Dissipation  Rate  Normalized  to 
((VAR)/C)3/2  Versus  Radar  Range,  Turbulence 
Outer  Scale  Length  are  X  -  0,  5,  4  km.  Radar 
full  half-power  beamwidtn  is  1°  and  pulse  volume 
depth  is  200  m.  Ratio  of  pulse  volume  to  depth  at 
top 
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If  >ni\  sni  i  trum  variance  data  is  available,  tiieri  ubscrvaUnn  •  >f  the  variation  •>( 
Doppler  spectrum  variance  with  range  (as  in  Figure  :>b)  may  identify  a  reasonable 
outer  scale,  thus  allowing  for  estimation  of  the  eddy  dissipation  rate.  In  vigorous 
storms,  such  as  thunderstorms,  such  a  technique  may  not  be  feasible.  Further¬ 
more,  because  of  the  storm  structure,  power  spectrum  analysis  of  the  Doppler 
mean  for  the  purpose  of  estimating  the  turbulence  outer  scale  may  be  of  little  use. 
For  such  cases,  it  may  be  possible  to  establish  a  generic  turbulence  outer  scale 
(sort  of  a  mean  storm  value)  perhaps  distinguished  by  storm  structure,  maximum 
equivalent  refleet  lvity  factor,  or  other  observable  storm  feature,  which  will  allow 
for  important  conclusions  concerning  turbulence  intensity  to  be  drawn.  1’his  con¬ 
cept  is  supported  by  Doviak  et  nl.  1  who  indicate  that  the  variance  of  the  Doppler 
spectrum  does  not  increase  dramatically  with  range  for  storm  ranges  of  50  to 
120  km.  1'his  suggests,  after  comparison  with  Figures  na  and  5b,  that  turbulence 
may  frequently  be  the  major  contribution  to  Doppler  spectrum  variance,  and  that 
on  the  average  the  turbulence  outer  scale  may  be  1  to  3  km.  This  result  would  be 
in  agreement  with  general  aircraft  observations. 

Consider  the  following  exaggerated  example.  An  outer  scale  length  of  2  km  is 
assumed  in  an  environment  where  X  actually  varies  from  1  to  4  km.  The  Doppler 
spectrum  variance  is  allowed  the  two  extreme  values  of  1  and  3fi  m".'ser~.  'Iable  1 
outlines  the  estimated  eddy  dissipation  rate  values  (e*'5),  a  quantity  typically 
used  as  a  turbulence  severity  indicator.  In  each  case  the  severity  estimate  derived 
from  the  generic  turbulence  outer  scale  value  estimates  w'ell  the  two  true  possible 
turbulence  intensities.  Thus,  use  of  a  reasonable  estimate  of  turbulence  outer  scale 
may  allow  for  accurate  classification  of  turbulence  severity  when  only  Doppler  spec¬ 
trum  variance  information  is  available.  For  ranges  greater  than  about  20  km, 
such  a  method  would  always  be  superior  to  the  commonly  used  relation  of  Frisch 
and  Clifford,  10  which  as  observed  in  Figure  fia,  will  always  underestimate  the 
eddy  dissipation  rate. 

Finally,  it  must  be  realized  that  all  the  results  derived  here  have  been  for  a 
particular  hypothetical  radar  having  a  pulse  volume  depth  of  200  m  and  full  half¬ 
power  beamwidth  of  1°.  Variations  in  pulse  volume  depth  should  noticeably  modify 
the  results  in  only  those  situations  where  pulse  volume  depth  is  the  largest  pulse 
volume  dimension.  Figure  7  displays  the  ratio  ^jqq/^Oo'  of  normalized  eddy 
dissipation  rate  from  two  radars  having  the  same  1°  beamwidth,  but  differing 
pulse  volume  lengths  of  100  and  200  m.  Perfect  tracers  were  assumed;  however, 
the  conclusions  are  applicable  to  all  precipitation  environments  considered  here. 

Note  that  Ejqq  =  e/((VAR)  /C)^^  and  for  X^  =  l  km  would  simply  be  the  X  =  l  km 

11.  Doviak,  R.J.,  Sirmans,  D. ,  Zrnic,  D. ,  and  Walker,  C.  R.  (1978)  Considera¬ 
tions  for  pulse  Doppler  radar  observations  of  severe  thunderstorms, 

J.  Appl.  Meteor.  17:189-205. 
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curve  in  Figure  fia.  Thus  e  ^qq  's  the  equivalent  curve  for  a  radar  having  a  100  m 
pulse  volume  depth.  A  ratio  of  1  at  some  range  means  that  the  two  curves  would  be 
colocated  at  that  range  and  both  radars  would  measure  the  same  Doppler  spectrum 
variance  and  eddy  dissipation  rate.  A  ratio  of  1  for  all  ranges  means  the  curves 
for  the  two  radars  would  be  identical.  The  plot  indicates  significant  difference 
only  for  ranges  less  than  about  10  km.  T  ius,  the  curves  presented  in  Figures  5 
and  (>  should  be  applicable  to  1°  beamwidth  radars  at  ranges  greater  than  about 
20  km  and  having  typical  pulse  volume  depths  of  100  to  300  m.  Additionally,  a 
change  of  10  to  20  percent  in  beamwidth  would  probably  result  in  a  similar  change 

1/3 

in  measured  Doppler  spectrum  variance.  The  resulting  change  in  e  (propor- 
1  /  2 

tional  to  (VAR)  '  )  should  be  much  less.  Thus  the  results  presented  here  should 
be  applicable  to  most  meteorological  radars  of  similar  dimensions. 

3.  AIRCRAFT  AND  RADAR  DATA  ANALYSES 

The  material  presented  so  far  deals  with  the  contribution  to  Doppler  spectrum 
variance  from  precipitation  in  a  given  turbulent  environment.  This  variance  has 
been  shown  to  be  directly  dependent  upon  the  two  basic  turbulence  parameters.  The 
eddy  dissipation  rate  is  the  quantity  of  interest  since  it  describes  the  intensity  of 
turbulent  air  motions  in  the  wind  field.  However,  Figures  a  and  9  demonstrate 
that  some  knowledge  of  the  turbulence  outer  scale  is  required  before  we  can 
successfully  estimate  this  quantity.  It  was  stated  that  turbulence  in  thunderstorms 
may  be  localized;  there  are  a  few  patches  wherein  the  turbulence  is  roughly  homo¬ 
geneous.  The  following  material  shows  results  of  first  analyses  of  aircraft 
thunderstorm  gust  velocity  data  which  appear  to  support  ’his  concept  and  yield 

estimates  of  both  e  and  A  . 

<  ) 

The  follow  :.g  data  were  acquired  on  May  2d,  10?fi  by  an  Air  Force  F-4  instru¬ 
mented  aircraft.  Two  penetrations,  separated  by  approximately  30  minutes,  were 
made  at  a  height  of  about  10,  000  ft  ( 4 8 7 f 1 .  3  m)  through  a  portion  of  a  storm  com¬ 
plex  located  40  km  south  of  Norman,  Oklahoma.  Aircraft  measurements  were 
taken  every  0.  1  second.  The  first  run  (run  1)  was  from  210°.  The  second  run 
(run  2)  passed  through  roughly  the  same  storm  region,  but  from  80°  and  30  minutes 
later.  Aircraft  speed  was  at  198  m  sec  during  both  runs.  The  discussion  will 
center  around  the  analysis  of  the  vertical  gust  velocity  data. 

Figures  8  (a  through  e)  and  9  (a  through  e)  are  plots  of  the  vertical  gust  veloc¬ 
ity,  derived  gust  velocity,  environmental  temperature,  aircraft  normal  accelera¬ 
tion,  and  barometric  altitude  during  the  two  runs.  The  prominent  feature  in  run  1 
is  the  apparent  strong  upward-moving  current  between  52.310  and  52540.  It  is 
about  <)  km  wide  and  is  characterized  by  a  warmer  temperature  than  the  environment 
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Figure  8.  Aircraft  Time  Series  Data  for  Early  Aircraft  Run 
of  (a)  Vertical  Gust  Velocity,  (b)  Derived  Gust  Velocity, 

(c)  Environmental  Temperature,  (d)  Normal  Acceleration, 

(e)  Barometric  Altitude.  Analyzed  segments  are  labeled  1 
through  7 


during  entry  and  increased  normal  acceleration  of  the  aircraft.  Also  note  the  large 
spikes  in  the  vertical  gust  velocity  of  a  somewhat  periodic  nature,  suggesting  the 
possible  presence  of  waves.  Similarly,  run  2  exhibits  two  apparent  upward-moving 
air  currents  centered  at  53775  and  53910  sec  respectively.  The  period  53555 
through  53645  exhibits  strong  periodic  fluctuations.  Increased  aircraft  normal 
acceleration  and  spikes  in  vertical  gust  velocity  are  also  noticed  here. 

Casual  observation  of  these  two  plots  suggests  that  each  gust  velocity  time 
sequence  may  be  segmented  into  a  number  of  discrete  local  turbulence  zones.  In 
run  1,  we  focus  on  three  zones  labeled  3,  5,  and  7  in  Figure  8,  part  (a),  corre¬ 
sponding  to  the  visually  estimated  turbulence  classes  moderate,  heavy  within 
updraft  region,  and  light,  respectively.  Similarly,  in  run  2  (Figure  9,  part 
(a)),  the  time  series  is  segmented  into  zones  1,  5,  6,  9,  7,  and  8  respectively. 
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Figure  9.  Aircraft  Time  Series  Data  for  Late  Aircraft  Run 
of  (a)  Vertical  Gus£  Velocity,  (b)  Derived  Gust  Velocity, 
(c)  Environmental  Temperature,  (d)  Normal  Acceleration, 
(e)  Barometric  Altitude.  Analyzed  segments  are  labeled  1 
through  9 
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Kaoh  gust  velocity  tune  scrii'S  «  a»  l<  d  to  power  spei  tr:d  analysis  (PSA). 

The  analyses  were  conducted  on  segments  of  the  original  (lata  sets.  Segments 
consisted  of  one  or  more  of  the  individual  /ones  and  are  labeled  1  through  7 
(run  1,  Figure  8)  and  1  through  l1  (run  2,  Figure  ft)  respectively.  Kach  segment 
was  linearly  detrended  before  being  subjected  to  spec  t  rum  analysis.  Although 
detrending  forces  the  power  spectral  density  to  go  to  zero  at  zero  frequency,  the 
spectral  outer  scales  are  generally  welt  beyond  the  influence  of  this  scale  removal 
action,  and  should  represent  a  true  change  of  spectral  slope.  Figures  10  and  11 
show  the  power  spectral  density  (PSD)  plotted  as  a  function  of  radian  frequency. 

An  approximate  space  settle  is  also  included.  Observations  front  run  1  will  be 
presented  first. 

It  is  observed  that  the  spectra  may  be  culssified  according  to  where  the  sharp 
changes  in  slope  (knee)  of  the  curves  occur.  Spectra  c,  d,  and  f  exhibit  apparent 
slope  changes  at  scales  less  than  1  km.  whereas  spectra  a,  b,  e,  and  g  appear  to 
have  the  knee  located  at  scales  greater  than  2  km.  It  is  believed  that  segment  7 
(spectrum  g)  would  also  have  shown  a  small  outer  scale  length  if  the  time  record 
52550  through  52555  had  not  been  included.  Reference  to  Figure  8  part  (a)  shows 
that  spectra  a,  b,  and  e  are  strongly  biased  by  inclusion  of  the  upeurrent  region  as 
a  turbulent  fluctuation  rather  than  as  a  true  feature  of  the  larger  storm  wind  struc¬ 
ture.  The  updraft  region  thus  acts  as  a  large,  very  energetic  eddy  which  stronglv 
2 

biases  the  strength  (V  )  of  the  turbulent  field  and  introduces  extra  energy  into  the 
long  wavelength  region  of  the  PSD  curves.  In  the  sense  that  the  upeurrent  does  not 
qualify  as  a  truly  turbulent  feature,  it  is  inappropriate  to  blithely  include  it  in 
PSA.  Spectra  d  and  f  show  the  power  spectra  derived  from  the  fluctuations  imbed¬ 
ded  in  the  positive  and  negative  slope  portions  of  this  upeurrent  region  It  is  seen 
that  these  spectra  are  similar  to  spectrum  c.  Referring  back  to  Figure  4b,  it  is 
observed  that  the  knee  of  the  spectral  curve  well  represented  the  turbulence  outer 
scale  length  under  the  ideal  conditions  of  homogeneous  isotropic  turbulence.  It  is 
felt  that  here  the  knee  may  be  a  reasonable  measure  of  this  parameter. 

Now  consider  run  2.  These  spectra  are  not  as  easily  categorized.  It  is  clear 
that  segment  4,  which  includes  all  data  as  turbulence,  exhibits  the  largest  appar¬ 
ent  outer  scale  value  and  generates  large  spectral  densities  at  long  wavelengths. 
Again,  this  results  from  inclusion  of  the  updraft  features  as  turbulent  gusts.  The 
lower  power  but  prominent  peak  near  .08  km  in  spectrum  a  results  from  a  0.  5  sec 
ripple  riding  on  the  gust  velocity  trace.  These  fluctuations  are  too  fine  to  be  seou 
in  these  highly  reduced  traces.  In  spectrum  f  a  long  (approximately  40  sec)  mod¬ 
ulation  is  the  likely  cause  of  the  spectral  peak  near  4.  G  km.  It  is  instructive  to 
compare  spectra  a,  c,  e,  and  f.  Note  that  spectra  a,  e,  and  f  have  a  well  defined 
knee.  However,  when  these  three  regions  are  combined  into  one  (3),  the  resulting 
spectrum  does  not  exhibit  this  clear  slope  change  at  short  scales.  In  this  instance. 
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2  3 

Figure  10.  Turbulence  Power  Density  (m  /sec  )  Versus  Frequency  (rad  sec) 
for  Segments  1  through  7  tn  Figure  8.  Part  (a)  corresponds  to  segment  1, 

(b)  to  segment  2,  etc.  .Apparent  space  scale  at  top 
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POWER  SPECTRUM  DENS  POWER  SPECTRUM  DENS. 


segment  f>  appears  to  aet  similarly  to  ail  uprurrent  region  for  the  eomposite  seg¬ 
ment.  The  apparent  turbulence  outer  scales  for  run  2  for  the  individual  segments 
are  in  the  range  .4  »«>  2.  4  km.  I  heSe  data  indicate  that  inclusion  of  storm  struc¬ 
ture,  notiturbulent  fluctuations,  and  eoinposites  of  local  turbulent  patches  <  an 
strongly  influence  the  result  ing  apparent  turbulent  power  spectra  derived  from  the 
data. 

Table  2  provides  the  mean  square  velocity,  estimated  outer  smile  length,  and 

estimated  eddy  dissipation  rate  for  these  various  spectra.  Note  that  the  outer 

scale  lengths  are  considerably  larger  for  run  2  than  for  run  1.  This  difference 

may  be  a  result  of  sampling  nonident  iral  storm  regions,  but  may  also  indicate 

a  temporal  evolution.  The  almost  synergistic  effect  of  combining  local  turbulence 

zones  into  single  larger  units  is  dramatically  demonstrated  bv  'lie  air  motion  var- 

2 

innce  V“.  Surprisingly,  the  eddy  dissipation  rates  for  the  individual  zones  of  the 

two  spectral  sets  are  very  similar.  Within  notiupdraf:  regions,  both  runs  show  e 

2  -3 

in  the  range  13  to  3 f  10  cm  see  .  The  largest  t  are  found  within  the  updraft  regions 

2  -3 

and  range  from  490  to  180,5  cm  See  .  These  values  indicate  that  nonupdraft  re¬ 
gions  would  be  characterized  as  moderate  in  strength,  while  the  updraft  regions 

,  12 

are  heavy  to  extreme. 


Table  2.  Kstimated  Turbulence  Parameters 


•Sped  rum 
Segment 

\*o. 

i 

2 

3 

4 

5 

fi 

7 

8 

2.  2.  2. 

V  (m  /s  1 

12.  19 

17.  73 

0.77 

7.  of. 

9.  25 

2.  81 

0.  41 

(Early  Period) 

A.  (km) 

o 

0.  I>2 

0.  82 

0.  95 

0.  5 

4 

0.  38 

,  2,  3* 

€(c;n  /sec  ) 

2749 

4823 

29 

18  85 

282 

490 

AT2/  2 /  2A 

V  (m  /sec  ) 

0.  534 

4.  15 

3.34 

32.  38 

5.  84 

1.  33 

2.39 

12.  38  (Late 

Period) 

A  (km) 

o 

1.  18 

1.  8 

2.48 

2 

1.  13 

0.  41 

2.  08 

2  3 

e(cm  /sec  ) 

13 

188 

99 

283 

54 

380 

847 

Finally,  it  is  instructive  to  correlate  the  aircraft  and  radar  data  to  determine 
if  the  radar  is  properly  measuring  the  turbulent  contribution  to  Doppler  spectrum 
variance.  This  must  be  ensured  before  radar  variance  data  are  interpreted  in 
terms  of  turbulent  field  parameters.  Before  presenting  the  results  in  detail. 
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however,  ills'  radar  data  and  the  technique*  applied  to  the  r:.d«r  an  !  aircraft  data 
to  prepare  them  for  the  rorrela'  inti  an  ilysis  will  be  described. 

The  radar  data  were  interpolated  •  -sito  a  (  a'  tcsiati  (x,  v,  '/)  grid,  "t  he  hori¬ 
zontal  surfaces  were  separated  vertically  by  Az  0.  f.  km.  The  grid  point  separa¬ 
tion  within  the  surfaces  was  also  Ax  Ay  0.  5  km.  A  horizontal  surface  near 
4  km  agl  represents  the  surface  through  which  the  aircraft  flew  and  radar  data 
fields  for  this  level  will  be  shown.  The  Doppler  spectrum  variances  were  obtained 
by  use  of  the  Fast  Fourier  Transform  on  or  iginal  complex  time  series  radar  data. 
The  Doppler  spectra  were  then  analyzed  using  an  objective  thresholding  scheme 
which  automatically  determined  the  spectrum  noise  level,  isolating  the  true  mete¬ 
orological  signal  from  which  the  spectral  mean  and  variance  were  determined. 

Figures  12  and  13  show  the  equivalent  reflectivity  factor  Z  (dBZ)  for  the  early 
and  late  aircraft  runs,  respectively.  The  storm  was  propagating  towards  the  N'E 
in  both  cases.  The  portions  of  the  aircraft  track  that  were  eventually  found  to  be 
best  correlated  with  the  radar  data  are  shown.  In  the  early  run  (Figure  12)  the 
aircraft  flew  through  one  cell  which  had  a  Z(  >  3d  dBZ,  and  was  part  of  a  storm 
complex.  The  actual  length  of  aircraft  track  also  observed  by  radar  is  about  11  km 
in  length  and  is  shown  by  the  line  segment  in  this  plot.  The  late  period  (Figure  13) 
shows  the  aircraft  flew  through  a  storm  consisting  of  one  large  storm  cell  of 
Z  >  3G  dBZ.  The  corresponding  aircraft  track  supported  by  radar  data  is 
about  38  km  long.  The  maximum  Zp  (<40  dBZ)  values  in  both  storms  are  not 
indicative  of  severe  weather;  however,  the  cells  are  strong  enough  to  require  nor¬ 
mal  avoidance  by  aircraft.  Figures  8,  part  (a)  and  9,  part  (a)  indicate  that  the 
aircraft  flew  through  at  least  one  updraft  region  in  each  run.  With  only  one  radar 
scanning  nearly  horizontally,  actual  vertical  velocities  cannot  be  measured.  How¬ 
ever,  information  from  plots  of  radar  radial  velocity,  shown  in  Figures  14  and  15, 
when  combined  with  the  Zp  plots  do  suggest  the  updraft  positions. 

Since  the  storms  were  moving  NE,  the  Zp  cells  would  normally  be  expected  to 
lie  just  upwind  (SW)  of  their  associated  updrafts.  Updrafts  in  horizontal  radial 
velocity  plots  generally  appear  as  regions  having  a  magnitude  near  zero.  Both 
radial  velocity  plots  show  regions  of  marked  decrease  of  radial  velocity  magnitude 
in  the  vicinity  of  the  high  cells  penetrated  by  the  aircraft.  In  the  late  period  the 
aircraft  penetrated  the  southern  end  of  a  velocity  minimum,  which  is  suggestive  of 
an  updraft  region.  In  the  early  period,  the  aircraft  penetrated  a  broad  region  of 
marked  decreased  magnitude,  but  not  a  well  defined  minimum.  These  regions 
agree  well  with  the  aircraft  updraft  positioning  (shown  as  U  in  Figures  14  and  15) 
determined  by  the  aircraft  track. 

Because  we  are  primarily  interested  in  the  turbulent  contribution  to  Doppler 
spectrum  variance  and  wish  to  correlate  the  radar  and  aircraft  Doppler  spectrum 
variance  estimates,  other  factors  which  broaden  the  radar  variance  estimates 
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must  be  removed.  At  the  low  elevation  angles  scanned  (<8°),  the  only  significant 
contributors  other  than  turbulence  are  the  transverse  shears  of  the  radial  wind. 
These  were  estimated  over-  the  level  of  aircraft  flight  and  are  shown  in  Figures  Id 
and  17,  and  Figures  18  and  19,  for  the  early  and  late  periods  respectively. 

The  vertical  shear  of  the  radial  wind  for  the  early  period  shows  maximum 

-3  -1 

shear  values  of  about  7  X  10  sec  .  Curiously,  the  only  Zg  cell  which  does  not 
show  a  shear  maximum  in  close  proximity  is  the  cell  penetrated  by  the  aircraft. 
Because  of  the  spectral  processing  techniques  applied,  large  transverse  shear 
values  found  outside  vigorous  storm  regions  may  generally  be  considered  represen¬ 
tative  of  the  environmental  wind  and  not  simply  a  result  of  low  spectrum  power. 

The  horizontal  shear  of  radial  velocity  plot  (Figure  17)  shows  a  markedly  periodic 
pattern  of  positive  and  negative  shear  maxima  through  the  entire  region  of  moderate 
storm  intensity  (Z^  >  24  dBZ).  The  structure  has  a  spacing  of  about  3  km  per 
maximum.  Considering  each  adjacent  pair  as  a  single  unit  (for  example,  each  pair 
could  represent  a  circulation,  or  vortex  eddy)  suggests  that  the  storm  complex  has 
a  basic  small  scale  structure  of  about  0  km.  This  scale  is  much  larger  than  the 
0.  5  km  grid  spacing  and  probably  represents  the  dominant  structure  scale  in  the 
storms. 

For  the  late  period  the  vertical  shear  also  has  a  maximum  value  of  about 
-3  -1 

7  X  10  sec  .  Also,  although  the  most  significant  shear  is  colocated  with  the 
single  large  Zg  cell,  strong  vertical  shear  is  again  observed  outside  the  storm 
cell.  The  horizontal  shear  shows  the  pattern  of  positive  and  negative  maxima, 
although  not  as  vividly  as  in  the  early  period.  The  small  scale  storm  structure 
here  is  closer  to  7  km  in  size.  The  shear  estimates  displayed  in  Figures  1G  through 
19  were  used  to  compute  the  shear  contribution  to  Doppler  spectrum  variance. 

These  contributions  were  then  subtracted  from  the  Doppler  spectrum  variance 
values,  resulting  in  variance  estimates  which  were  dominated  by  turbulent  con¬ 
tributions. 

The  turbulence  produced  Doppler  spectrum  variance  estimates  are  shown  in 

Figures  20  and  21.  The  early  period  (Figure  20)  shows  that  the  variance  is  gen- 

2-2  ? 
erally  <6  m  sec  outside  the  storm  (Z  <  20  dBZ)  and  usually  closer  to  2  m 

-2  *  o  -9 

sec  .  The  high  Z  cells  within  the  storm  complex  exhibit  variance  >4  m  sec" 
e  2-2 

and  are  generally  closer  to  8  m  sec  .  Within  the  storm  complex  in  general 

2  -2 

(20  dBZ  <  Z  <36  dBZ)  the  variance  values  are  in  the  range  2  to  6  m  sec  .  The 

®  2-2 
late  period  (Figure  21)  similarly  shows  spectrum  variance  to  be  <4  m  sec  outside 

the  storm  (<20  dBZ),  >8  sec  ^  within  the  high  Zg  cell,  and  2  to  6  sec  ^  else¬ 
where  in  the  storm. 

As  with  any  general  rule  there  are  exceptions.  Good  examples  where  large 
Doppler  spectral  variance  exists  outside  high  Zg  cells  are  seen  at  (-4,  35)  in 
Figure  20,  and  (-22,  33)  and  (-2,  47)  in  Figure  21,  Noting  that  large  variance 
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Turbulence  Induced  Doppler  Spectrum  Variance  (m  /sec  )  on  Horizontal 
'ing  Early  Aircraft  Run 


generally  represents  energeiie  turbulent  conditions  (Figure  (1),  we  see  that  these 
failures  show  that  hazardous  zones  may  frequently  be  found  outside  high  Zp  cells 
and  outside  what  one  would  generally  consider  the  outer  storm  boundary.  Next  we 
return  to  the  task  of  comparing  these  radar  spectrum  variance  estimates  with  the 
equivalent  aircraft  derived  values. 

The  aircraft  gust  velocity  data  consist  of  a  time  series  of  vertical  gust  com¬ 
ponents  separated  in  distance  by  about  0.02  km.  To  obtain  equivalent  Doppler 
spectrum  variance  estimates  from  aircraft  data  the  following  procedure  was  fol¬ 
lowed.  First,  a  continuous  time  series  of  gust  velocity  components  was  formed  by 
taking  the  individual  detrended  segments  (3,4,  (>,  and  7  for  the  early  run  [Figure  8]; 
and  1,  a,  (i,  9,  7  and  8  for  the  late  run  [Figure  9]  )  and  coupling  them  to  form  a  single 
time  series  of  "true”  gust  velocity.  These  series  are  shown  in  Figures  22  and  23. 
Second,  a  Gaussian  filter  representing  the  radar  beam  pattern  function  was  applied 
to  each  time  series  to  obtain  a  time  series  of  equivalent  Doppler  spectrum  variance 
for  each  run.  These  variance  estimates  are  thus  estimates  of  what  the  radar  would 
have  measured  if  it  ooserved  these  aircraft  vertical  gust  velocity  components  in  a 
field  of  uniform  reflectivity.  Next,  to  correlate  this  data  set  with  the  actual  radar 
variance  data  shown  in  Figures  20  and  21,  the  aircraft  data  were  interpolated  to  a 
set  of  grid  points  along  the  aircraft  track.  The  grid  points  were  separated  by 
0.5  km.  This  technique  results  in  aircraft  and  radar  data  which  have  undergone 
similar  beam  filtering  and  interpolation  to  Cartesian  grid  locations. 

Before  proceeding  further,  however,  we  must  take  note  of  deficiencies  in  this 
method.  First,  the  aircraft  equivalent  variance  is  derived  from  the  vertical  turbu¬ 
lent  gust  velocity  while  the  radar  spectrum  variance  is  derived  from  essentially  the 
horizontal  component  of  turbulent  gust  velocity.  Second,  for  the  aircraf.  measure¬ 
ment  the  radar  beam  pattern  function  was  applied  to  a  line  of  aircraft  data,  but  to 
a  volume  defined  by  the  radar  pulse  volume  for  the  radar  measurement. 


Figure  22.  Time  Series  of  "True"  Vertical  Gust  Velocity  for  Early 
Aircraft  Run 
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Figure  23.  Time  Series  of  "True1'  Vertical  Gust  Velocity  for  Late  Aircraft  Run 


Also  note  that  use  of  a  0.  5  km  grid  point  system  introduces  error  since  radar  data 
can  be  interpolated  to  a  grid  point  from  regions  through  which  the  aircraft  did  not 
fly.  One  is  forced  to  use  some  sort  of  grid  system,  however,  since  the  aircraft 
position  actually  fluctuates  off  the  mean  flight  path  by  hundreds  of  meters.  Last, 
because  the  radar  was  not  slaved  to  the  aircraft  location,  temporal  evolution  of 
turbulent  fields  between  the  times  when  the  radar  and  aircraft  sampled  the  same 
general  region  could  have  occurred.  Thus,  differences  between  the  radar  and  air¬ 
craft  Doppler  spectral  estimates  are  not  unexpected. 

The  best  correlated  Doppler  spectrum  variance  data  sets  for  aircraft  and 
radar  are  shown  in  Figures  24  and  25.  The  early  data  set  (Figure  24)  has  a  corre¬ 
lation  coefficient  of  0.891  and  the  late  period  (Figure  25)  has  a  correlation  coef¬ 
ficient  of  0.  821.  Considering  the  deficiencies  of  the  technique  mentioned  above, 
the  degree  of  correlation  is  quite  remarkable  and  suggests  that  the  sampled  turbu¬ 
lent  fields  may  be  nearly  isotropic.  In  both  plots  the  region  of  largest  spectrum 
variance  is  in  best  agreement,  suggesting  that  conditions  approaching  isotropicity 
are  more  likely  in  regions  of  strong  turbulence  than  in  regions  of  light  turbulence. 

The  fact  that  these  data  correlate  well  suggest  that  turbulence  parameters  determined 
from  the  radar  data  alone  may  generally  estimate  well  the  actual  environmental 

values.  Last,  it  is  noteworthy  that  the  minimum  radar  variance  is  generally  about 

2-2  2-2 
2  m  sec  while  the  corresponding  aircraft  values  are  closer  to  0.25  m  sec 

This  may  result  from  deficiencies  in  the  technique  mentioned  above,  or  indicate 


that  additional  noise  or  small  scale  storm  structure  may  still  be  contaminating  the 
radar  data. 

It  is  unclear  whether  the  discrepancies  in  these  variance  plots  result  from 
deficiencies  in  the  technique  as  noted  above,  or  from  an  incomplete  understanding 
of  how  the  turbulent  precipitation  motions  are  mapped  into  Doppler  spectrum  vari¬ 
ance.  It  is  also  unclear  whether  the  agreement  seen  here  would  be  retained  on  a 
radar  pulse  volume  scale,  rather  than  the  0.  5  km  grid  scale  employed.  To  answer 


these  questions,  more  complete  data  sets  where  the  aircraft  measures  all  three 
gust  velocity  components  and  where  the  radar  pulse  volume  is  slaved  to  the  air¬ 
craft  location  are  required. 

The  results  discussed  here  will  allow  for  increased  accuracy  in  Doppler  radar 
estimation  of  turbulence  intensity,  and  will  have  practical  application.  As  stated 
earlier,  the  two  main  contributions  to  Doppler  spectrum  variance  are  the  trans¬ 
verse  shear  of  the  radial  wind  and  turbulence.  The  Doppler  radar  can  estimate  the 
shear  component  by  observing  the  variation  in  azimuth  and  elevation  of  the  mean 
Doppler  velocity  between  successive  radials.  Thus,  the  turbulent  contribution  can 
be  estimated.  With  some  knowledge  of  the  turbulence  outer  scale  and  storm  reflec¬ 
tivity  factor,  a  reasonable  estimate  of  the  eddy  dissipation  rate  may  be  obtained. 
This  information  can  then  be  combined  with  the  known  response  characteristics  of 
aircraft  for  identification  of  regions  of  hazard  to  aircraft.  Use  of  these  results 
with  incoherent  radar  is  significantly  more  difficult,  however.  The  incoherent 
radar  cannot  estimate  the  Doppler  mean  velocity,  and  therefore  will  not  allow  for 
the  determination  of  the  transverse  shear.  Thus,  the  turbulent  variance  contribu¬ 
tion  may  not  be  well  known. 

13 

Atlas  and  Srivastava  proposed  a  method  for  estimating  an  atmospheric  struc¬ 
ture  function  with  incoherent  radar  by  forming  combinations  of  signal  returns  from 
successive  pulse  volumes  along  a  given  radial.  They  showed  that  the  mean  squared 
difference  of  Doppler  velocity  between  the  adjacent  pulse  volumes  can  be  related  to 
the  Doppler  spectral  variances  by 

(V2  '  Vl)2  -  T  «VAR)12(pi  +  P2)  -  <VAR>lPl  -  (VARN2P2) 

(33) 

where  X  is  the  radar  wavelength,  (VAR)^  and  (VAR)2  are  the  Doppler  spectral  var¬ 
iances  (sec  )  from  pulse  volume  numbers  1  and  2  having  total  spectrum  powers 
Pj  and  P2,  and  (VAR)^2  is  the  variance  of  a  composite  Doppler  spectrum  obtained 
from  adding  the  signals  from  the  individual  pulse  volumes.  Unfortunately,  this 

relation  produces  a  highly  smoothed  estimate  of  the  true  atmospheric  structure 

14 

function,  requiring  very  careful  interpretation;  it  is  also  related  to  the  environ¬ 
mental  wind  shear. 


13.  Atlas,  D. ,  and  Srivastava,  R.C.  (1971)  A  method  for  radar  turbulence 

detection.  IEEE  Trans,  on  Aerospace  and  Electronics,  AES-7,  No.  1, 
pp  179-187. 

14.  Sychra,  J.  (1974)  Fluctuation  Spectra  and  Velocity  Structure  Function. 

Tech.  Rpt.  33,  Lab.  Atmos.  Prob. ,  U.  of  Chicago. 


A  measured  Doppler  radial  estimate  may  be  considered  as  a  sum  of  a  large- 
scale  background  wind  (V  ),  turbulent  component  (Vt),  and  an  additional  random 
error  and  noise  bias  (e),  as 

V  =  Vw  +  V(  +  e  .  (34) 

With  this  convention  the  left  side  of  Eq.  (33)  becomes 

(V2  -  Vl)2  =  (V2w  -  Vlw)2  +  (V2t  -  Vlt>2  +  <e2  '  el)2  (35) 

where  the  second  right  hand  term  is  actually  the  true  Doppler  estimate  of  the  pulse 
volume  filtered  turbulence  structure  function.  If  two  azimuthal  pairs  of  pulse  vol¬ 
umes  ([Vg,  Vj]  and  iVg,  Vgl)  are  combined  in  this  manner,  then  in  a  region  where 
the  mean  turbulent  air  velocity  over  the  region  occupied  by  the  pulse  volumes  varies 
little, 

(V3  -  V,)2  -  (V2  -  V,)2  .  (V3w  -  Vlw)2  -  (V2.  -  Vlw)2  .  (36) 

Using  Eq.  (36),  a  linear  variation  in  large-scale  wind  can  be  fitted  to  the  data, 
thereby  providing  an  estimate  of  the  large-scale  shear.  If  a  quadratic  variation  in 
the  wind  field  is  desired,  two  such  sets  01  Eq.  (3  >)  may  be  used. 

Obviously  this  technique  is  prone  to  error  where  -no  gradient  of  mean  turbulent 
velocity  is  large.  However,  at  the  ranges  of  interest  (perhaps  60  to  200  km)  most 
of  the  turbulent  energy  is  mapped  into  the  Doppler  spectrum  variance.  For  this 
case,  the  radar  estimate  of  the  turbulent  structure  function  is  small,  and  its  varia¬ 
tion  is  not  significant,  allowing  Eq.  (36)  to  be  useful.  The  validity  of  this  relation 
is  independent  of  the  actual  magnitude  of  the  air  turbulent  intensity.  Therefore, 
the  above  relations  may  be  applicable  within  local  patched  of  turbulence  and  may 
return  reasonable  estimates  of  transverse  radial  wind  shear,  thus  permitting  an 
estimate  of  eddy  dissipation  rate  to  be  obtained.  Usefulness  of  this  technique  can 
only  be  established  after  comparison  of  radial  shear  estimates  from  Doppler  and 
incoherent  techniques  using  real  data. 


4.  CONCLUSIONS 

The  results  of  this  theoretical  investigation  demonstrate  that  useful  and  under¬ 
standable  relationships  exist  between  the  mean  and  variance  of  Doppler  spectra  and 
turbulence  in  precipitation  environments.  It  is  found  that  Doppler  spectrum  vari¬ 
ance  and  estimated  eddy  dissipation  rate  are  strongly  dependent  upon  the  precipitation 


58 


n 


environment  for  ranges  less  than  about  20  km,  and  also  for  rases  where  'he  turbu¬ 
lence  outer  scale  length  is  less  than  about  0.  5  km.  The  Doppler  spectrum  variance 
and  eddy  dissipation  rate  are  essentially  independent  of  the  turbulence  outer  scale 
length  when  the  maximum  pulse  volume  dimension  is  less  than  about  one-half  the 
outer  scale  value.  As  range  increases,  both  parameters  are  strongly  dependent 
upon  range  and  outer  scale  length.  Once  the  largest  pulse  volume  dimension 
exceeds  the  outer  scale  length,  however,  Doppler  spectrum  variance  and  eddy 
dissipation  rate  are  independent  of  range  and  dependent  only  upon  the  turbulence 
outer  scale.  Aircraft  data  analyses  suggest  that  aircraft  gust  time  series  data 
should  be  viewed  as  a  composite  of  segments  of  local  turbulent  patches  wherein 
the  turbulence  parameters  may  be  essentially  constant.  Agreement  between  air¬ 
craft  and  radar  data  support  the  theoretical  concepts,  but  point  to  the  need  for 
more  definitive  data  sets.  While  use  of  equivalent  reflectivity  factor  may  be  suf¬ 
ficient  to  classify  the  precipitation  environment,  further  investigation  to  determine 
reliable  methods  for  remotely  estimating  the  turbulence  outer  scale  is  needed  if 
reasonable  estimation  of  the  turbulence  eddy  dissipation  rate  is  to  be  accomplished. 
Nonetheless,  good  estimates  of  the  severity  of  turbulence  may  be  obtained  by  use 
of  a  reasonably  guess  of  the  turbulence  outer  scale  length,  and  should  prove  useful 
in  detecting  hazardous  zones  within  storms. 
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